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ABSTRACT 

We study the total mass-density profile for a sample of 14 fast-rotator early-type galaxies (stellar 
masses 10.2 < log M*/Mq < 11.7). We combine observations from the SLUGGS and ATLAS 3D sur¬ 
veys to map out the stellar kinematics in two-dimensions, out to a median radius for the sample of four 
half-light radii R e (or 10 kpc), and a maximum radius of 2.0-6.2 R e (or 4-21 kpc). We use axisym- 
rnetric dynamical models based on the Jeans equations, which allow for a spatially varying anisotropy, 
and employ quite general profiles for the dark halos, and in particular do not place any restriction 
on the profile slope. This is made possible by the availability of spatially extended two-dimensional 
kinematics. We find that our relatively simple models provide a remarkably good description of the 
observed kinematics. The resulting total density profiles are well described by a nearly-isothermal 
power law Ptot( r ) °c r -7 from i? e /10 to at least 4 R e , the largest average deviation being 11%. The 
average logarithmic slope is ( 7 ) = 2.19 ± 0.03 with observed rms scatter of just tr 7 = 0.11. This 
scatter out to large radii, where dark matter dominates, is as small as previously reported by lensing 
studies around r « R e / 2, where the stars dominate. Our bulge-halo conspiracy places much tighter 
constraints on galaxy formation models. It illustrates the power of two-dimensional stellar kinematics 
observations at large radii. It would now be important to test the generality of our results for different 
galaxy types and larger samples. 

Keywords: galaxies: elliptical and lenticular, cD - galaxies: formation — galaxies: kinematics and 
dynamics — galaxies: structure 


1. INTRODUCTION 

One of the pillars on which our entire paradigm of 
galaxy formation rests is the fact that dark matter of an 
unknown nature dominates the mass budget of the Uni¬ 
verse (Blumenthal et al. 1984). Without dark matter, 
the primordial gas would not be able to collapse suffi¬ 
ciently quickly within the center of dark matter halos to 
form the galaxies we observe (White & Rees 1978). 

Total mass-density profiles of spiral galaxies were ob¬ 
tained decades ago, given the simple geometry of their 
spiral disks and the availability of ionized (Rubin et al. 
1980) and neutral gas (Bosma 1978), which provides a 
kinematical tracer easy to measure and model (see re¬ 
view by Courteau et al. 2014). The measurements in¬ 
dicated flat circular rotation curves and consequently 
nearly-isothermal p tot oc r -2 profiles. 

The situation is very different for early-type galaxies 
(ETGs: ellipticals and lenticulars), which by definition 
lack well-defined spiral disks and are cold-gas poor. For 
ETGs one usually has to rely on expensive observations 
of the stellar kinematics and more complex dynamical 
models. This has restricted most studies to radii not 
much larger than the half-light radius R e . 

A general consensus has emerged for the mass distri¬ 
bution of the ETGs inner parts. Long-slit observations 


of two different samples of « 20 ETGs revealed rota¬ 
tion curves to be nearly flat with nearly-isothermal mass 
profiles, as in spiral galaxies, within the median radius 
r ss 2 R 0 sampled by the kinematics (Gerhard et al. 2001; 
Thomas et al. 2011). This agrees with strong gravita¬ 
tional lensing studies finding nearly-isothermal slopes, 
with small scatter, for the total galaxy density profile 
of 73 ETGs, at a typical radius of r ss R e /2 (Auger et al. 
2010). These central slopes are similar to those mea¬ 
sured for group/cluster-scale profiles using X-ray mod¬ 
eling (e.g. Humphrey & Buote 2010) or stacked weak- 
lensing technique (e.g. Gavazzi et al. 2007). 

The largest detailed study of dark matter in galaxy 
centers was based on the integral-field observations of 
the volume-limited ATLAS 3D sample of 260 ETGs (Cap¬ 
pellari et al. 2011). It inferred a median dark matter 
fraction of just /dm {Re) = 13%, within a sphere of ra¬ 
dius R e , over the full sample (Cappellari et al. 2013a). 
This shows that studies restricted to the central regions 
of ETGs mainly measure the stellar mass distribution. 

To explore the region where dark matter dominates, 
one needs to reach at least ~ 4 R e . Existing studies 
at these radii targeted one galaxy at a time. They 
used observations of extended HI disks (Weijmans et al. 
2008), deep stellar kinematics at a few sparse locations 
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Table 1 

Sample of fast rotator early-type galaxies and measured parameters 


Name 

(u 

D 

(Mpc) 

(2) 

log R c 
(") 
(3) 

ttmax 

(kpc) 

(4) 

-Rmax/ Re 

(5) 

<7>r 

(0.1-li?e) 

(6) 

<7>r 

(l-4i?e) 

(U 

<7>r 

(0.1-4Re) 

(8) 

(|A|)r 

(0.1-4i? e ) 

(9) 

MGE ref. 

(10) 

A jam/Xloess 

(11) 

mediandiSVimsI) 

(%) 

(12) 

NGC 0821 

23.4 

1.60 

13.2 

2.92 

-2.16 

-2.36 

-2.23 

0.016 

C06 

1.09 

8.5 

NGC 1023 

11.1 

1.68 

11.3 

4.39 

-2.20 

-2.19 

-2.20 

0.009 

S09 

0.98 

12.3 

NGC 2768 

21.8 

1.80 

19.7 

2.95 

-2.10 

-1.82 

-2.01 

0.025 

C06 

1.02 

5.1 

NGC 2974 

20.9 

1.58 

10.8 

2.80 

-2.22 

-2.49 

-2.30 

0.015 

C06 

1.10 

14.9 

NGC 3115 

9.5 

1.54 

10.1 

6.19 

-2.28 

-2.53 

-2.37 

0.028 

E99 

1.10 

10.7 

NGC 3377 

10.9 

1.55 

8.7 

4.64 

-2.22 

-1.75 

-2.05 

0.046 

C06 

0.99 

8.4 

NGC 4111 

14.6 

1.08 

3.8 

4.46 

-2.12 

-2.16 

-2.13 

0.043 

— 

0.85 

7.7 

NGC 4278 

15.6 

1.50 

9.0 

3.76 

-2.19 

-2.45 

-2.29 

0.028 

C06 

1.09 

7.8 

NGC 4473 

15.3 

1.43 

8.4 

4.18 

-2.12 

-2.29 

-2.18 

0.015 

C06 

1.20 

6.4 

NGC 4494 

16.6 

1.69 

8.0 

2.04 

-2.18 

-2.52 

-2.26 

0.045 

S13 

1.06 

9.2 

NGC 4526 

16.4 

1.65 

18.2 

5.13 

-2.21 

-2.29 

-2.24 

0.014 

C06 

0.97 

9.0 

NGC 4649 

17.3 

1.82 

21.4 

3.86 

-2.10 

-2.34 

-2.19 

0.020 

S13 

1.05 

8.9 

NGC 4697 

11.4 

1.79 

12.5 

3.66 

-2.16 

-2.34 

-2.23 

0.023 

S13 

0.99 

7.1 

NGC 7457 

12.9 

1.56 

6.9 

3.05 

-1.82 

-2.23 

-1.96 

0.036 

C06 

1.33 

6.1 


Note. — Column (1): galaxy name; Columns (2)-(3): distance and half-light radius from Cappellari et al. (2011); Columns (4)-(5): 
Maximum radius _R ma x sampled by the stellar kinematics; Column (6): Average logarithmic slope ( 7 ) r = A log ptot /A log r of the total 
mass profile in the interval 0.1R e < r < R e ; Column (7): as in Col. 6 for R e < r < min(4R e , Rmax); Column (8): as in Col. 6 for 
0.1.Re < r < min(4i? e , i?max); Column (9): Average absolute deviation between ptot and the best-fitting power-law (in dex); Column (10): 
Reference for the MGE model: C06=Cappellari et al. (2006), S09=Scott et al. (2009), S13=Scott et al. (2013), E99=Emsellem et al. (1999) 
Column (11): Quality of fit, where XJAM is measured from the JAM models and XLOESS from the smoothed data. In both cases the 
reference is the original SLUGGS data, excluding values with zero final weight. This ratio approximates x 2 /DOF, but is insensitive to 
the normalization of the kinematic uncertainties. Column (12): Median of the absolute relative deviations |JVrms| = IVrms/K^2 ESS — 1| 
between the original and the smoothed SLUGGS Ums- 


(Weijmans et al. 2009; Forestell & Gebhardt 2010; Mur¬ 
phy et al. 2011), globular cluster kinematics (Napolitano 
et al. 2014) and planetary nebulae (Romanowsky et al. 
2003; de Lorenzi et al. 2008, 2009; Das et al. 2011; Napoli¬ 
tano et al. 2011; Morganti et al. 2013). 

The fact that galaxies were modeled using different 
techniques or kinematic tracers, and for a mix of ETGs 
with axisymmetric and triaxial shapes, resulted in a still 
unclear picture of the global trends. Tentative conclu¬ 
sions suggest a trend with massive ETGs having nearly 
isothermal total mass profiles and flat circular velocities 
(see review by Gerhard 2013), while lower mass ones have 
more slowly rising dark matter profiles and correspond¬ 
ingly falling circular velocity profiles (also Romanowsky 
et al. 2003). 

The situation is changing with the availability of ex¬ 
tended stellar kinematics for significant samples of ETGs 
(Brodie et al. 2014; Raskutti et al. 2014). Here we present 
the first detailed and fully homogeneous analysis using 
larges-scale stellar dynamics of a statistically significant 
sample of 14 ETGs (but see Deason et al. 2012 for less 
detailed models). The key novelty of this work is the 
availability of two-dimensional stellar kinematics out to 
a median radius for the sample of r « 4 R c from the 
SLUGGS survey (Brodie et al. 2014), which we com¬ 
bine with integral-field kinematics in the central regions 
(~ 1 R e ) from the ATLAS 30 survey (Cappellari et al. 
2011 ). 

2. SAMPLE AND DATA 
2 .1. Axisymmetric sample selection 

We study a subsample of the SLUGGS ETG sample, 
with kinematics from Arnold et al. (2014). We want 
a homogeneous set of galaxies and we use an axisym¬ 
metric dynamical modelling method. For this reason we 
only select ETGs classified as fast rotator in Emsellcm 
et al. (2011). We add the fast rotator NGC 3115, without 
ATLAS 30 kinematics. Our sample of 14 galaxies is given 
in Table 1 (stellar masses 10.2 < log M*/M 0 < 11.7 from 


Cappellari et al. 2013a). 

Fast rotators were shown by Krajnovic et al. (2011) to 
have kinematical axes, within about 1 R e , extremely well 
aligned with the photometric major axes measured at the 
much larger radii ~ 3 R e sampled by the SLUGGS kine¬ 
matics. This alignment, for the statistically-significant, 
ATLAS 30 sample, shows that fast rotator ETGs, as a 
class, must be axisymmetric out to the region sampled 
by our models. Significant deviations from axisymme- 
try would produce a broad distribution of kinematical 
misalignment that is clearly ruled out by the ATLAS 30 
data. The SLUGGS kinematics of fast rotators is also 
generally consistent with axisymmetry and shows small 
kinematical misalignment. 

Some genuine deviations from a bi-symmetric velocity 
field do exist in the SLUGGS data. Notable examples 
are NGC 4494 in the maps of Arnold et al. (2014) and 
NGC 4473 in Foster et al. (2013). However, these devi¬ 
ations are more likely due to unrelaxed substructure in 
the stellar halo, than to non-axisymmetry in the relaxed 
stellar distribution. In fact, the tight alignment between 
ATLAS 30 stellar kinematics and photometry, as well as 
the complete lack of photometric twist, is extremely un¬ 
likely in a triaxial configuration (e.g. Binney 1985). In 
this work we attribute deviations from the axisymmetric 
assumption to either sub-structure or systematic prob¬ 
lems in the data, and simply try to remove their effects 
from our models. 

2 .2. Photometry and mass models 

We use the Multi-Gaussian Expansion (MGE) (Em- 
sellern et al. 1994; Cappellari 2002) to parametrize both 
stellar and dark halo distributions. The MGEs of 8 galax¬ 
ies come from Cappellari et al. (2006) and one from Scott 
et al. (2009). They include deep /-band photometry out 
to 5-10/? e - Three SDSS r-band MGEs come from Scott 
et al. (2013). The U-band MGE for NGC 3115 comes 
from Emsellcm et al. (1999). The MGE for NGC4111 
was redone on the i-band SDSS image using the Python 
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Figure 1. Symmetrization and cleaning of SLUGGS data. Left panels show the stellar Vims = VV 2 + <r 2 , while right ones show the 
mean stellar velocity V. (a) original linearly-interpolated data; (b) bi-symmetrized data; (c) bi-symmetrized and LOESS-smoothed data. 
Note the good prediction, in the central parts, of the actual SAURGN kinematics in Figure 2; (d) best-fitting JAM model. The right-panel 
shows the JAM prediction for the mean velocity using the same model parameters as in the left-panel, adopting for all Gaussians a^ = an 
and k = 1.20. We tried different a^/an ratios 3 but found best fits for a^ & an (see Cappellari 2008 for definitions). 


version 1 of the MGE-fitting method of Cappellari (2002), 
to avoid the r-band saturation (see Scott et al. 2013, sup¬ 
plementary material). 

2.3. Stellar kinematics 

The kinematics 2 3 for the inner (~ 1 R e ) region of our 
galaxies was observed with the SAURON integral-field 
spectrograph and was homogeneously extracted for the 
ATLAS 30 survey as described in Cappellari et al. (2011). 
For 10/14 galaxies the observations, taken as part of 
the SAURON survey, were already presented in Em- 
sellem et al. (2004). The kinematics for the outer regions 
were observed with the Keck/DEIMOS multi-slit spec¬ 
trograph by the SLUGGS survey. They were presented 
and described in Arnold et al. (2014). 

Briefly, in both cases the kinematics were extracted in 
pixel space using the pPXF method (Cappellari & Em- 
sellern 2004). The ATLAS 30 data covered the H/? and 
Mg b spectral region (4800-5380 A) and were extracted 
using the full MILES stellar library (Sanchez-Blazquez 
et al. 2006) as templates. While the SLUGGS data cov¬ 
ered the Ca II triplet spectral region (8480-8750 A) and 
were extracted using the stellar library of Cenarro et al. 
( 2001 ). 

3. METHODS 

3.1. Symmetrization and outliers removal 

A key step of this analysis is the removal of potential 
outliers from the SLUGGS kinematics. These may have 

1 Available from http://purl.org/cappellari/software 

2 Available from http://purl.org/atlas3d 

3 This is done by replacing T> in equation (38) of Cappellari 
(2008) with [D + (1 — c^b^q^a^./R 2 ] (Note 9 of the arXiv version 
of that paper). 


significant effects on the modeling results, owing to the 
relatively sparse sampling of the kinematics at large radii. 

When fitting kinematics to models, one can enforce 
symmetry during the kinematic extraction (e.g. Cap- 
pcllari & Emsellem 2004). We adopt here a conceptu¬ 
ally similar idea, which is easier to apply to irregularly 
sampled or discrete kinematics (see also Morganti et al. 
2013). 

At every observed location ( Xj , yj ) on the sky, where 
the x-axis coincides with the kinematic major axis PA^m 
(taken from Krajnovic et al. 2011), we calculate the rms 
velocity U r 2 nsj - = V) 2 + er 2 from the SLUGGS data, where 
Vj is the mean stellar velocity and cr ? is the velocity dis¬ 
persion. We then generate a new set of bi-symmetric 
positions by replicating four times the observed values 
as follows: 

Xj —/ (Xj,—Xj,Xj,—Xj) 

yj (yj,yj,-yj,-yj) (i) 

Kms,j ^ (^rms,j 5 Kmsj 1 Kms ,j ■> Kms ,j ) 

The measurement uncertainties are increased by 2x. 
Given the bi-symmetry of the model, the y 2 statistics, 
describing the quality of the fit, is unchanged. Equa¬ 
tion (1) applies to the even velocity moments like U rm s 
and a. While for the odd velocity moments like V, the 
last line is replaced by: 

( 2 ) 

The resulting kinematics measurements are then 
smoothed with the (Python version) of the two- 
dimensional LOESS algorithm of Cleveland (1979), 
implemented 1 by Cappellari et al. (2013b) (see example 
in Figure 1). 

Comparison between the smoothed/cleaned V(n2 ESS 
and the original one suggests the SLUGGS error being 
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Figure 2. Dynamical models from two-dimensional stellar kinematics. In each panel, the top panel shows the observed stellar V^ms, while 
the bottom panel shows the JAM model prediction. The spatial scale of the top panels is in units of the half-light radius R e . Contours of 
the MGE surface brightness in 1 mag steps are overlaid. The data include DEIMOS/SLUGGS stellar kinematics at large radii. They were 
bi-symmetrized and LOESS smoothed as described in the text (see Figure 1). The locations sampled by the data are shown as black dots. 
The SAURON/ATLAS 3D data, used at smaller radii, are enclosed by the magenta lines. 
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Figure 3. Same as in Figure 2 for the second half of the galaxy sample. 
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underestimated due to systematics, with typical errors 
of ~ 10% (Table 1). The original errors would indicate 
non-smooth kinematics, which seems inconsistent with 
the smooth galaxy images. 

In the following models we do not fit the original kine¬ 
matics, but the symmetrized/smoothed one. This has 
the added benefit of making the agreement between data 
and model easy to visually assess. We need to verify that 
our models capture the global features of the data, and 
that our results are not driven by a few deviant values. 
This is important in situations where data systematics 
may be present and relying entirely on y 2 statistics may 
be misleading. 


constrained but are not of interest here. They are just 
‘nuisance parameters’, marginalized over to derive the 
total mass profiles studied here. The parameters are: (i) 
the inclination i; (ii) the anisotropy j3 z = 1 — cr 2 /cr^, 
with a z and ctr the stellar dispersion in cylindrical coor¬ 
dinates, for the MGE Gaussians with 07 < f? e ; (iii) the 
anisotropy for the remaining Gaussians at larger radii; 
(iv) the stellar (M/L) stars ; (v) the break radius of the 
dark halo, constrained to be 10 < r s < 50 kpc; (vi) the 
halo density p s at r s ; and (vii) the dark halo slope a for 
r <C r s . 

4. RESULTS 


3.2. Weighting and matching of the two datasets 

Another issue for the modeling is the fact that the 
ATLAS 30 observations consist of many more data points 
with smaller uncertainties, which completely dominate 
the x ' 2 estimate. However here we want our dark 
halo determinations to be especially constrained by the 
SLUGGS data, which sample the region where the halo 
dominates. Similarly to Morganti et al. (2013), we in¬ 
crease the size of the ATLAS 30 kinematic uncertainties 
so that, for a good fit, the two datasets provide an equal 
contribution to the y 2 . We leave the SLUGGS uncertain¬ 
ties unchanged, to retain properly normalized confidence 
levels for our model parameters. 

The SLUGGS U rms at the SAURON locations tend to 
be lower than the measured SAURON values. We find a 
median offset of 11%, which is larger than the 5% level 
we consider unavoidable between independent datasets. 
This offset was noted by Arnold et al. (2014), but its 
source is unknown. The ATLAS 30 data agree on aver¬ 
age with hundreds of independent determinations from 
the literature (fig. 8 of Cappellari et al. 2013a), suggest¬ 
ing the SLUGGS data may be offset with respect to the 
optical literature. Here we simply multiply the SLUGGS 
U rms to fit, for each galaxy, the interpolated SAURON 
data at the same locations. This is the standard kine¬ 
matics we fit with our models. However, importantly, 
we have also run all our models with the SLUGGS data 
alone, and confirmed that the slopes of the total mass 
profiles agree with those of our standard models. 

3.3. Dynamical models 

We model the ATLAS 30 and SLUGGS stellar kine¬ 
matics using the (Python version) of the axisymmetric 
Jeans anisotropic modeling (JAM) method 1 (Cappellari 
2008). The approach is the same used in Cappellari 
et al. (2013a) for the ATLAS 30 data alone. It employs 
a Bayesian method with constant (i.e., ‘ignorant’) priors 
on all parameters. 

The key difference between this work and previous stel¬ 
lar dynamical studies of dark halos in ETGs is the fact 
that we place virtually no constraint on the halo profile 
parameters. The halo is assumed to be spherical, but is 
described by a generalized Navarro et al. (1996) profile 
(gNFW) with free normalization, inner slope and break 
radius: 

Our models have 7 free parameters. Some are poorly 


4.1. The models describe the data well 

The first result is the simple fact that the models pro¬ 
vide a good description of the stellar kinematics of all 
the modeled galaxies over the full field (Figure 2, 3 and 
Table 1). This was not expected. It would have been 
natural if we had employed e.g. the more general orbit- 
superposition method (e.g. Cappellari et al. 2006), which 
is fully described by thousands of parameters. However 
our models have just six non-linear parameters and one 
scaling factor ( M/L) stais . 

Moreover the fits look similarly good even assuming 
a power-law halo and a constant-anisotropy stellar body 
(four nonlinear parameters and one scaling factor). The 
fact that four parameters are able to consistently de¬ 
scribe all features of the two-dimensional maps, for the 
full set of 14 galaxies, then indicates that (i) the (cleaned) 
SLUGGS data are reliable and (ii) the model assump¬ 
tions provide a good description of the dynamics and 
mass distribution of the real galaxies. 

An alternative interpretation for the good fits would 
be that the anisotropy and dark matter variations, and 
data systematics, conspire to mimic the simple orbital 
structure and mass distribution assumed by the mod¬ 
els. This would be a realistic possibility when studying 
a single galaxy, given the non full generality of the JAM 
models, but such a conspiracy is unlikely for such a large 
and diverse set of galaxies. 


4.2. Isothermal profiles with small scatter to Ji? e 

The second and main result of this work is that all 
14 modeled fast-rotator ETGs have a nearly-isothermal 
Ptot ex r~ 2 total density distribution from R e /10 out to 
the median radius of 4i? e sampled by this study (Fig¬ 
ure 4(a)). The total mass-density profiles 4 , marginalized 
over all nuisance parameters, are tightly constrained by 
the data. In the whole range 0.1i? e < r < 4 R e the 
profiles are well described by a power-law p tot oc r -7 
with the largest average deviation of 11%. The corre¬ 
sponding average logarithmic slope is ( 7 ) = 2.19 ± 0.03 
for the sample, with an rms scatter of just <r 7 = 0 . 11 . 
We do not detect any significant correlation of the slope 
with R e , stellar mass, or stellar velocity dispersion. For 
0.1i? e < r < R e the average slope is (7) = 2.15 ± 0.03 
with a 1 = 0.10, while for R e < r < 4R e we find 
(7) = 2.27 ± 0.06 with ct 7 = 0.23. 


1 Computed from the axisymmetric MGEs as 
M i exp[-r 2 /(2af)] erf [r^/l-q? / ( qj c 
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Figure 4. Profiles of the total mass-density distribution, (a) Measured profiles for the 14 modeled galaxies with SLUGGS+ATLAS 30 
data. The profile for each galaxy was plotted for 100 realizations randomly drawn from the posterior distribution of the model parameters, 
to illustrate the random model uncertainties. Three lines with p oc r (NFW inner slope), p oc r~ 2 (isothermal) and p oc r ~ 3 are also 
shown, (b) Cosmologically-motivated profiles (these are models E from Cappellari et al. 2013a). These were computed by attaching, to the 
stellar density of the ATLAS 30 galaxies, spherical NFW dark halos with masses predicted by the abundance matching technique in such a 
way that the models fit the stellar kinematics. These models naturally predict a nearly isothermal total-mass profile out to r > 10 R e . (c) 
For comparison the purely-stellar profiles of the same ATLAS 30 galaxies in (b) are shown, (d) Comparison between our density profiles 
and published ones. The profiles with outline enclose the allowed range of published profiles, while the colored bands are realizations from 
our models posterior. 


Our inner-profile slope and scatter are in excellent 
agreement with those ( 7 ) = 2.08 ± 0.03 with <r 7 = 0.16 
found around r ss R e /2 using strong lensing (Auger et al. 
2010 ). 

Figure 4(b) shows that the observed trend is consistent 
with what one would predict for the whole ATLAS 3D 
sample, for cosmologically-motivated uncontracted NFW 
halos. In Figure 4(c), the stellar profiles are very different 
from the total ones at the radii we sample: they have 
slopes pstars oc r ~ 2 around r ~ R e /2 (fig. 2 of Cappellari 
et al. 2013a) but fall off more steeply than p s t ars oc r -3 
around r ss 4i? e - 

Figure 4(d) compares our total profiles with published 
ones for NGC 0821, NGC 2974, NGC 4494, NGC 4649, 
and NGC 4697 (from Forestell & Gebhardt 2010; Weij- 
mans et al. 2008; Morganti et al. 2013; Das et al. 2011; 
de Lorenzi et al. 2008, respectively). The ptot{r) was de¬ 
rived from the circular velocities v c [r ), assuming spher¬ 
ical symmetry. In 4/5 cases the agreement is excellent, 
with our statistical uncertainties enclosing the published 
profiles over the full radial range. Notable is the tight 
agreement for NGC 2974, where the v c (r) was directly 
measured from a regular HI disk. We believe the dis¬ 


agreement for NGC 4494 may be due to the inclusion in 
Morganti et al. (2013) models of the strong asymmetry 
in the SLUGGS data at r > 100". We excluded those 
data from our fits, but including them would improve the 
agreement. 

5. CONCLUSIONS 

We combine the integral-held stellar kinematics from 
the ATLAS 30 survey, within ~ 1 R e , with the two- 
dimensional stellar kinematics from the SLUGGS survey, 
out to a median radius of about 4i? e , and a maximum 
radius of 2.0-6.2 R e , for a sample of fast rotator early- 
type galaxies consistent with axisymmetry. We construct 
the first statistically significant set of detailed axisym- 
rnetric dynamical models of the two-dimensional stellar 
kinematics out to those large radii, where dark matter 
dominates. 

We find that the galaxies’ dynamics is well represented 
by few relatively simple assumptions. The models tightly 
constrain the total density profiles, which closely approx¬ 
imates the isothermal form p t ot oc r^ 2 , from R e /10 out 
to the median radius of 4i? e sampled by the data, with 
remarkably little scatter. The observed total mass distri¬ 
bution is not a generic prediction of ACDM and provides 
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constraints on the models (e.g. Remus et al. 2013; Dutton 
& Treu 2014). 

Our sample highlights the importance of similar stud¬ 
ies on larger samples of galaxies, to provide a much 
needed benchmark for galaxy formation models. For this, 
studies like the present one, using DEIMOS on Keck or 
MUSE on the VLT, can be complemented with mod¬ 
els of shallower data, but for much larger samples like 
MaNGA (Bundy et al. 2015). To be most useful, sam¬ 
ples need to include both spiral and early-type galaxies. 
These should be modeled in a fully homogeneous way by 
describing the kinematics of both the stars and the gas 
within the same gravitational potential. Large radio sur¬ 
veys of the neutral HI gas like Apertif (Verheijen et al. 
2008) can complement MaNGA at the largest radii. The 
future availability of large kinematic datasets indicates a 
bright future for the systematic investigation of the mass 
profiles in galaxies. 
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